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Recently a formulation of overlap fermions at finite density based on an analytic continuation 
of the sign function was proposed. We study this proposal by analyzing the energy and number 
densities for free fermions as a function of the chemical potential and the temperature. Our results 
show that overlap fermions with chemical potential give rise to the correct continuum behavior. 
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1. Introduction 

In recent years finite temperature and finite density physics has seen a lot of attention from the 
lattice QCD community. Although so far mainly traditional fermion formulations such as staggered 
or Wilson fermions are used, it is an important issue to get also exactly chiral fermions, in particular 
the overlap operator |jlj|, ready for QCD thermodynamics. 

In |Q] Bloch and Wettig have suggested a formulation of overlap fermions which includes 
a chemical potential. The proposal is based on introducing the chemical potential as the 4-th 
component of the gauge field [[|] and an analytic continuation of the sign function into the complex 
plane, which requires new techniques for the evaluation of the overlap operator [|J]. 

In a recent paper [||] we have studied thermodynamical properties of the new proposal by 
analyzing the energy density for the free case as a function of the chemical potential and the tem- 
perature. In this contribution we augment this analysis by computing also the number density 
and comparing it to the known result for a free fermion gas. Our study indicates that the pL- and 
/"-dependence for both quantities, energy and number density, approach the continuum behavior 
correctly. 

2. Setup of the calculation 

The overlap Dirac operator D{pL) for fermions with a chemical potential pL is given as [Q] 

D(ji) = i[l-y 5 sign// Gu)], 

H(jl) = y$[l -aDw(jl)], (2.1) 

where the sign function may be defined with the spectral theorem. Dw(jti) denotes the Wilson 
Dirac operator with chemical potential, 

D w (n)x, y = + — }s x>y (2.2) 
la c?4 J 

3 r 



For later use we distinguish between the lattice spacing a for the spatial directions (= 1,2,3 di- 
rections) and the temporal lattice constant an,. We use periodic boundary conditions in the spatial 
directions and anti-periodic boundary conditions for time. The chemical potential \i is introduced 
in the standard exponential form 

For jit = the Wilson Dirac operator is 75-hermitian, i.e., /s/VW/s = Av(0)^> implying that 
H(0) is a hermitian matrix. When the chemical potential /I is turned on, 75-hermiticity no longer 
holds, and H(/j.) is a non-hermitian, general matrix. This fact implies: First, the eigenvalues of 
H(/j.) are no longer real and the sign function for a complex number has to be defined in the spectral 
representation of sign//(/x). Second, the spectral representation has to be formulated using left- 
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and right-eigenvectors. For the sign function of a complex number we use the analytic continuation 
proposed in [||] and define the sign function through the sign of the real part 



sign (x + iy) = sign (x) . 
For analyzing thermodynamic properties we study the following observables: 
1. The energy density 

1 dlnZ 



(2.3) 



2. The number density 



V dp 



T d\nZ 
V d^i 



(2.4) 



(2.5) 



Here Jj? is the Hamiltonian of the system, denotes the number operator, and J3 = l/T is the 
inverse temperature (we set the Boltzmann constant equal to 1). The derivatives in the definition of 
e(ju) are taken such that fin = c = const. 

The continuum result for the subtracted thermodynamical quantities of free massless fermions 



can be calculated using the definition ( |2.4j ) and (|2.5|). The final expressions together with the 
logarithm of the partition function can be found in textbooks (see, e.g., [^]). 

In the path integral formalism the inverse temperature j3 is given by the extent in 4-direction, 
which on the lattice has the form /3 = N^a^. Thus the derivative d/dfi in (2.4) is expressed as 
N^ 1 d / da^. The partition function Z is given by the fermion determinant detD which we write as 
the product over all eigenvalues X n of D(pL). After some algebra one obtains 



e(jl) 



1 d In detD 



VN 4 da 4 
1 d In detD 

VN^an djj. 



£l4/i=C 



1 1 dX n 



a4jl=c 



1 y J_ dX n 

VN 4 a 4 V K dfl 



(2.6) 



(2.7) 



3. Evaluation of the eigenvalues 



For the evaluation of (2J5) and (2.7) we need the eigenvalues A of D(ju) in closed form. This 
is done in three steps: First, using Fourier transformation, we bring the Dirac operator for free 
fermions to 4 x 4 block-diagonal form, 



(3.1) 



with 



h 5 = l-Y j [l-cos(ap j )]-—[l-cos(a 4 (p4-in))] , 



"j 
h 4 



7=1 



a 4 



— sin(apj) for j = 1,2,3 , 

- — sin(a 4 (p 4 — in)) . 
a 4 



(3.2) 
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The spatial momenta are given by pj = Inkj/aN, where N is the number of lattice points in the 
spatial directions and kj = 0, 1 ...N — 1. The momenta in time-direction are p4 = 7l(2&4 + \)/a 4 N 4 , 
fct = 0,l...AT 4 -l. 

In the second step the spectral representation for the sign function of H is constructed. The 
left- and right-eigenvectors of H, L and rj are given by 



/ 



i -. zj 0) [// + sl] , l 2 = l K 2 ! [H + sl] , h = l^'[H-st] , U = l?>[H-st], 

.(o 

3 



f(°)l 



/(0)[ 



n = [H + st]^ , r 2 = [H+sl]r^ , r 3 = [H - st]^ , r 4 = [H-sl]rf ] 



(3.3) 



;(°) JO) 



The constant spinors Zj , r\ are (T denotes transposition) 



7 (0) 

(o) 



I 



(0)7 



(0) 
2 



c(l, 0,0,0) , I 
c (0,0, 1,0) , 4 0) 



r (0)r 

2 



c(0, 1,0,0), 
c (0,0, 0,1). 



(3.4) 



where c = (2^(5 + /js)) ! / 2 and^= y/j 2 + /jj. The eigenvectors obey Z,r ; - = The corresponding 
eigenvalues GCj of // are 

«i = «2 = +5 , Of3 = 0:4 = — s . (3.5) 
Using these eigenvalues and eigenvectors we can construct sign H using the spectral 
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sign// = £ sign (A/) r/Z; 



sign (5 J 



//. 



(3.6) 



Inserting this result in (2.1), the eigenvalues of the overlap operator at a given momentum are 
obtained as 



A± 



1 



sign (Jh 2 + hj)h 5 ± iVf? 



h 2 + h 2 5 



(3.7) 



where each of the two eigenvalues is twofold degenerate. In the spectral sums ( |2.6| ), ( |2.7D the label 
n runs over all momenta and the eigenvalues at fixed momentum as given in (3/7). The necessary 
derivative with respect to a 4 and pL is straightforward to compute in closed form, and the spectral 
sums (2.6), (2.7) can then be summed numerically. The argument of the sign function cannot 
become purely imaginary on a finite lattice, and no 5-like terms occur. We remark, that after taking 
the derivative with respect to a 4 , we set a = a 4 = 1, i.e., all the results are presented in lattice units. 



4. Results 

We begin the discussion of our results with lattices which have equal length in all four direc- 
tions, corresponding to zero temperature in the thermodynamic limit. We focus our attention on 
the subtracted energy density e(ju) — e(0) and the number density n p (pl) —n p (0). The results for 
various lattice sizes are shown in Fig. 1. According to [|6j], we expect the data (symbols in Fig. 1) 
to approach the continuum form 

eOx)-e(0) = i^ , MM)-» P (0) = At, (4.1) 
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Figure 1: The energy density g(/x) — e(0) and the number density n p (jx) — n p (0) as a function of jj, 4 and 
/I 3 , respectively. The symbols (connected to guide the eye) are for various lattice sizes, the dashed line is 
the continuum result. All quantities are given in lattice units. 

represented by dashed lines in Fig. 1. 

We see clearly that both quantities show almost linear dependence on \i 4 and /x 3 , respectively. 
The figure also demonstrates that with enlarging the lattice the finite size effects become smaller 
and the lattice results approach the continuum limit. The deviation from the linear dependence 
is due to the finite extension of the lattice. Furthermore, we can expect that we can see finite 
temperature effects too. 

To study the effects of temperature, we compare our data with the continuum results [Q] for 
the subtracted energy and the number density, now including the leading temperature corrections, 

u 4 T 2 u 2 u 3 uT 2 

e(H,T)-e(0,T) = ^2+^- , n p (n,T)-n p (0,T) = ^ + tt—. (4.2) 

Fig. 2 shows the lattice results for finite temperature, using data from 128 3 x Afy lattices with 
temperature T = N^ 1 (in lattice units). The approximate linearity of the zero-temperature case re- 
ceives corrections as the temperature is increased. The temperature-dependent terms have a square- 
and cubic -root behavior when we plot them as a function of pL A and n 3 , respectively. 

In order to study the temperature effect quantitatively, we compare the finite temperature re- 



sults to the continuum form (4.2). We have to take into account also finite size effects which we 
build in assuming contributions from higher terms of the expansion in }i. In the case of the energy 
density, we use the ansatz 

e( J u,r)-e(0,r) = C2M 2 + c 4 M 4 + c 6 M 6 + + ^(iO • (4-3) 

For the number operator we use an equivalent form with odd powers of jj,. The results of the fit for 
the data used in Fig. 2 and for the largest lattice of Fig. 1 are given in Table 1. 

As the table shows, the coefficients approach the continuum values when the temporal exten- 
sion is increased. In the case of the energy density, the discrepancy for N4 = 24 is approximately 
2% and 9% and for only 0.5% and 7%. For N4 smaller than 24, the discrepancy for both 
coefficients is larger, which is due to finite size effects. Nevertheless both, the }i. A - and /i 2 r 2 -terms, 
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Figure 2: The energy density e(ji, T) — e(Q,T) and the number density n p (jj., T) —n p (0, T) as a function of 
jX 4 and /J, 3 , respectively, for finite temperature lattices. All quantities are in lattice units. 



are well reproduced. The coefficients C(, and eg have magnitudes of <^(10~ 2 ), so their contribu- 
tion according to ( |4.3[ ) is irrelevant for 0.1 < pL (in lattice units). For the number density we find 
qualitatively the same behavior with only slightly different numerical values. 

Finally, we analyze the /j. -dependence for large values of the chemical potential, approaching 
the cutoff. In this case we compare our data with data obtained for the Wilson fermions and the 
continuum case. In Fig. 3 we plot the ratio (£(/x) — £(m)) /jU 4 - The figure shows that in a certain 
interval, from zero up to approximately 0.7 in lattice units, both lattice formulations give practically 
identical results and start to differ only beyond this already relatively large value of /J.. The second 
observation is that up to jJ. « 0.6 the lattice results differ very little from the continuum form and 
discretization errors grow large only for }i > 0.6. Also in a comparison with the data for fermions 
obtained by blocking [Q] one again finds good agreement. 



N 4 


N 4 2 /2 




c 4 


D(c 2 ) 


D(c 4 ) 


8 


0.007812 


0.010125 


0.03519 


29 % 


39 % 


12 


0.003472 


0.004125 


0.03178 


19 % 


25 % 


16 


0.001953 


0.002192 


0.02803 


12 % 


11 % 


24 


0.000868 


0.000947 


0.02587 


9 % 


2 % 


128 


0.000030 


0.000032 


0.02543 


7 % 


0.5 % 



Table 1: Results for fits according to (4.3). The spatial volume is always 128 3 . The temporal extension 
N4 is given in the first column. In the second column we list the corresponding value of N^ 2 /2 which is 
what one expects for the fit coefficient C2 in the third column. The coefficient C4 is expected to approach the 
constant value 1 / 4tt 2 = 0.02533. D{c2) and D(c4) are the relative deviations of C2 and 04 from the expected 
continuum values. 
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Figure 3: The ratio of (e(/x) — e(jtt)) / ji A as a function of jj.. The figure compares overlap fermions with 
Wilson results and the continuum. 

5. Summary 

We have analyzed thermodynamical quantities for the overlap operator at finite chemical po- 
tential obtained by analytic continuation of the sign function into the complex plane as suggested 
in The analysis was done for finite and zero temperature lattices. Fits of the data show that the 
expected continuum behavior is reliably approached and no trace of unphysical p} /a 2 terms was 
found. We conclude that the suggested introduction of chemical potential [Q] provides both, chiral 
symmetry and the correct description of fermions at finite density. 
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